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Pair  correlations  and  the  effects  of  dense  packing  are  included  to  extend  existing  results  on 
reflection  and  scattering  by  distributions  of  protuberances  (bosses)  on  rigid  or  free  base  planes. 
The  earlier  energy  conserving  forms  for  the  specularly  reflected  wave,  the  surface  impedance,  and 
the  differential  scattering  cross  section  per  unit  area,  are  obtained  in  terms  of  a  transform  of  the 
scattering  amplitude  of  an  isolated  boss.  Low-frequency  approximations  are  developed  with 
emphasis  on  the  roles  of  the  packing  density  and  multipole  coupling  effects,  and  explicit  results 
are  given  for  semi-elliptic  cylinders  and  hemi-ellipsoids.  For  lossless  bosses,  the  reflection 
coefficient  has  a  minimum,  and  the  incoherent  scattering  a  maximum,  at  the  packing  density 
corresponding  to  maximum  fluctuations  in  the  number  of  bosses  per  unit  area.  Multipole 
coupling  effects  may  be  misinterpreted  in  data  inversion  programs  as  changes  in  boss  shape;  if 
such  effects  are  not  included,  then,  e.g.,  hemispheres  may  be  mistaken  for  hemiellipsoids 
broadened  along  the  base  plane  and  shortened  along  the  normal.  - .  . 


PACS  numbers:  43.20.Fn,  43.20.Bi 


INTRODUCTION 

Various  physical  aspects  of  the  coherent  reflection  and 
incoherent  scattering  of  sound  by  rough  surfaces  are  exhibit¬ 
ed  by  random  distributions  of  protuberances  on  rigid  or  free 
base  planes.  For  sparse  uncorrelated  distributions  we  ob¬ 
tained1  energy -conserving  approximations,  general  enough 
to  account  for  the  dependence  on  angle  for  near-grazing  inci¬ 
dence  and  for  pseudo-Brewster  effects,  and  simple  enough  to 
analyze  special  shapes  1,2  and  to  reduce  data  in  experimental 
contexts.3 

The  representations'  of  the  reflection  coefficient  (R ) 
and  the  differential  scattering  cross  section  [Q )  in  terms  of  an 
impedance  (Z )  proportional  to  the  scattering  amplitude  ( /)  of 
one  protuberance  (boss)  in  isolation  and  to  the  number  ( p)  of 
bosses  per  unit  area,  are  restricted  to  sparse  packings.  At 
dense  packings  and  low  frequencies,  i.e.,  for  negligible  inco¬ 
herent  scattering,  we  supplemented  R  with  results  for  peri¬ 
odic  distributions4-3  for  the  case  of  one  propagating  mode 
(separation  of  bosses  small  compared  to  wavelength  A  ).  To 
investigate  domains  of  validity  for  the  sparse  and  periodic 
cases  for  parallel  cylindrical  protuberances  (analogs  of  the 
statistical  mechanics  one-dimensional  cases  of  the  sparse  gas 
and  crystal),  we  included  the  effects  of  pair  correlations.6 
Using  the  Zemike-Prins7  one-dimensional  liquid  state  pair 
function  pp(x)  for  minimum  separation  b  and  average  separ¬ 
ation  b  =p~'  of  centers,  plus  a  rapidly  converging  residue 
series  for  b<b  and  b~b,  led  to  forms  in  terms  of  pb  =  w  for 
the  full  range  0  <  w<  1,  from  sparse  gas  on  to  the  determinis¬ 
tic  periodic  limit.  The  present  paper  extends  the  earlier  work 
on  striated  surfaces6  to  facilitate  data  inversion  purposes, 
and  provides  analogous  results  for  rough  surfaces  of  arbi¬ 
trary  bounded  bosses  in  terms  of  a  pair  function />/>(R)  with  R 
as  the  separation  of  centers  of  pairs.  The  development  is 
based  on  recent  work  on  correlated  monolayers.8 
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We  obtain  the  original  form 1  for  A  (Z )  in  terms  of  a  more 
complete  impedance  proportional  to  a  transform  of /  that 
includes  packing  effects  and  correlations,  plus  a  correspond¬ 
ing  energy-conserving  approximation  for  Q.  We  construct 
explicit  low-frequency  results  (with  semi-ellipses,  and  hemi¬ 
ellipsoids  as  illustrations)  and  emphasize  the  roles  of  packing 
density,  and  multipole  coupling.  It  is  shown,  for  example, 
that  for  lossless  bosses,  the  specularly  reflected  energy  has  a 
minimum,  and  the  incoherent  scattering  a  maximum,  at  the 
packing  density  corresponding  to  maximum  fluctuations  in 
the  number  of  bosses  per  unit  area.  Another  application  re¬ 
lates  to  the  inversion  of  reflection  data  to  determine  boss 
shape:  e.g.,  because  of  multipole  coupling  effects  a  hemis¬ 
phere  corresponds  to  an  equivalent  hemi-ellipsoid  (broad¬ 
ened  along  the  base  plane  and  shortened  along  the  normal), 
and  such  effects  need  to  be  included  to  the  data  reduction 
procedure. 

For  hemispheres  with  symmetry  axis  perpendicular  to 
the  plane,  we  take  the  packing  density  ( w )  either  as  the  frac¬ 
tion  of  plane  covered  by  the  boss  bases  (of  radius  a),  or  more 
generally  as  the  fractional  area  covered  by  cocentered  larger 
acoustically  transparent  disks  (of  diameter  b>2a).  The  b 
disks  correspond  to  impenetrable  statistical  mechanica  par¬ 
ticles9"12  with  b  equal  to  the  minimum  separation  of  centers; 
equivalently,  b  is  the  radius  of  the  exclusion  region  contain¬ 
ing  the  center  of  only  one  base  of  radius  a.  The  additional 
parameter  b  for  the  minimum  separation,  or  the  correspond¬ 
ing  packing  density  w  =  ptrb  2/4,  enables  us  to  model  a 
broader  class  of  surfaces.  The  size  of  w,  from  near  zero  on  to 
wd  =;0.84  (an experimented  value  for  the  densest  packing  of 
identical  circular  disks)  determines  the  correlation  effects. 
More  generally,  the  base  of  the  boss  and  of  the  exclusion  disk 
are  not  taken  as  circular;  their  shapes  and  orientationm  the 
plane  are  specified  by  cocentered  curves  a(R)  and  b  (R)  that 
are  neither  similar  nor  similarly  aligned.  As  before,8  we  as¬ 
sume  that  b(R)  fully  determines  /?(R)  and  that  both  have  the 
same  reflection  and  inversion  symmetries  as  an  ellipse. 
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In  the  following  we  use  expressions  given  in  Refs.  1  and 
8  and  refer  to  earlier  equations  by  ( 1 :38),  etc.  For  brevity,  the 
recent  work  on  monolayers8  is  cited  for  many  detailed  re¬ 
sults  that  are  required.  We  work  with  forms  for  the  average 
multiple  scattering  amplitude  F  for  one  boss  in  terms  of  the 
isolated  value  /  and  in  terms  of  transforms  of /  (We  could 
construct  these  forms  by  superposing  results  for  the  symmet¬ 
rical  monolayer,  but  we  include  some  additional  relations 
F[f ]  to  facilitate  perturbation  extensions  to  more  general 
base  planes.)  For  detailed  applications,  we  consider  series 
decompositions  of  F  and  /  in  terms  of  the  same  scattering 
coefficients  A„  and  a„,  and  the  distribution  integrals 
(continuum  analogs  of  lattice  sums)  derived  for  the  symmet¬ 
rical  monolayer.8 

I.  COHERENT  REFLECTION 

As  before,8  we  use 

r  =  rr,  r  —  z2  +  R2,  f  =  z  cos  9  +  R(^)sin  9, 

R(<p)  =  x  cos  <p  +  y  sin  <p,  ( 1 ) 

as  well  as  the  direction  cosines  r  =  zy  4-  xa  4-  y  0  =  £,  z,  y, . 
Suppressing  e  we  work  with  a  pair  of  plane  waves, 
imaged  in  the  base  plane  z  —  0, 

<t>  =  e* r,  k  =  kk,  k  =  2ir/ A, 

k  =  r (90,  <plt)  =  f0, 

6  '  =  e,v  \  k'  =  kk', 

k'  =  Hit  -  9a,  <p0)  =  f £  =  f0. .  (2) 

We  write  k'  =  k  -  2kzz  =  k  -  2y£  with  y0  =  cos90 
=  Icos0o|.  For  cylinders^ with  generators  along  y,  we  set 
¥  =  ¥o<  0  =  7}  =  0,  and  R  =  x.  In  general  we  use  three-di¬ 
mensional  terminology;  if  two  forms  of  a  factor  arise,  we  list 
the  cylindrical  first. 

For  <b  incident  on  an  arbitrary  obstacle  at  the  phase 
origin,  we  write1,8  the  scattering  amplitude  asg(r,  k).  If  the 
obstacle  is  symmetrical  to  reflection  in  z  =  0,  then  g(r, 
k)  =  g( r',  k')  with  r'  =  r(jr  —  9,  <p).  The  scattering  amplitude 
for  the  corresponding  boss  on  a  rigid  (  +  )  or  free  (  —  )  base 
plane  excited  by  ±  <t>  ’  follows  by  superposition, 

/±  I?-  k')  =g{r,  k)  ±  git  k') 

=  git,  k')  ±  git,  k)  =/±  (?',  k).  (3) 

The  boss  amplitudes  ft  are  obtained  directly  from  known 
results  for  g,  as  discussed  originally  by  Rayleigh”  for  fine 
semicircular  cylinders.  Equivalently,  f+  or  /_  is  twice  the 
component  of  g  symmetrical  or  antisymmetrical  to  reflec¬ 
tion  in  z  =  0.  We  normalize  g\r,  k)  as  before, 1,8  so  that  for 
lossless  cases,  —  Re  g(k,  k)  equals  the  mean  value  of  |g(f,  k)|2 
over  all  directions  of  observation  r. 

We  restrict  consideration  to  the  half-space  z  >  0  (hence¬ 
forth  the  back  space).  The  choice  ±  <f> '  for  the  incident  wave 
for  the  corresponding  rigid  and  free  base  planes  allows  for  a 
single  representation  of  the  scattered  and  reflected  fields, 
and  for  the  general  suppression  of  subscripts  ±  .  Thus  the 
specular  scattering  theorem 14  for  either  case  has  the  form 


-  o()  Re/(k,  k ' )  =  <7a  +  a, ,  a,  =  1  /|r,  k)  | 2; 

_  4  Air 
a°~~k’  p-; 

«  pn/2  «  riT  pn/2 

= —  d9,  —  dtp  d9s\n9,  (4) 

7 1  J  -  IT/ 2  2  IT  J  —  TT  JO 

where  cra  and  at  are  the  absorption  and  scattering  cross  sec¬ 
tions;  ^|/|2,  the  mean  value  over  the  back  space  of  real 
directions  (0<<?<7r/2),  is  also  (by  symmetry)  the  mean  over¬ 
all  directions.  See  Ref.  1  for  additional  theorems  and  discus¬ 
sion  of  f±  for  arbitrary  bosses,  and  for  explicit  low-  and 
high-frequency  approximations  for  semicircular  cylinders 
and  hemispheres;  results  for  elliptic  semicylinders  are  given 
in  Ref.  2,  and  for  hemi-ellipsoids  in  Ref.  5  (p.  661,  IT). 

Similarly  for  </>  incident  on  a  configuration  of  obstacles 
(s  =  1,2,.. .A)  with  centers  at  R5  on  z  =  0,  we  specify  the  re¬ 
sponse  as  in  (8:7)  in  terms  of  the  multiple  scattering  ampli¬ 
tudes  Gs  determined  functionally  by  the  isolated  values  gs  as 
in  G  [g]  of  (8:8).  For  the  corresponding  boss  problem  we  work 
with  F  [/],  such  that  in  terms  of  G ’s  for  a  symmetrical  mono- 
layer,  F±  (r,  k')  =  G{ r,  k)  ±  G(r,  k')  =  G( r',  k')  ±  G  (r ',  k) 
equals  twice  the  symmetry  components  of  G.  Thus,  with 
e'V  R  Fs  as  the  amplitude  for  the  boss  at  Rs,  we  write  the 
solution  (excess  pressure)  for  the  configuration  of  bosses  as 

*±  =  ±<y  +  <t>  +  &  ±  =  ±<t>'  +  ¥£, 

^(r)  =  X  e'VB  f  RjF,(rr,  k'), 

S  Jc 

kc=krc;  rc  =  r(9c,  <pc),  (6) 

where  equals  (\/ir)fd9c  with  contour  as  for  H'q,  or 
(\/2ir)Sdfi(9c,  <pc)  with  contours  as  for  h Similarly, 

F, (?,£')=/, (f,  £') 

+  y2s'  J/,(f',  KtFslK,  kYK~k,R", 

R„=R, -Rs,  (7) 

w'here  1'  is  the  sum  over  ,s#t. 

For  a  homogeneous  ensemble  of  configurations  of  iden¬ 
tical  aligned  bosses,  we  proceed  as  for  (8:9)  to  obtain  the 
corresponding  ensemble  average  of  (6), 

( ^  >  =  p  j «*■  r{Fs{i' ,  k))s  J  d  Re'lk  -  “■ 1 R 

=  <b  2CF(k',  k); 

C  =p/ky0,  pir/k  2y0,  (8) 

where  F  =  (F, )  s  is  the  ensemble  averaged  amplitude  for  any 
fixed  boss.  Thus 

{¥)=  ±<y  +  <t>  +{%t)=  ±<y +{*£),  (9) 

with  the  coherent  reflected  wave  as  in  (1:37), 

<**>=*(1+2CF)  =  *I±!-. 

F=F(k,k')  =  F(k',k),  (10) 

in  terms  of  the  nomalized  impedance  Z.  The  same  form  also 
arises  for  periodic  surfaces  with  p~ 1  as  the  grating  spacing 
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for  parallel  cylinders4  or  as  the  area  of  unit  cell  for  the  rec¬ 
tangular  lattice5  and  spacings  small  compared  to  A  (one 
propagating  mode). 

To  obtain  F  =  (F, ) , ,  we  consider  the  ensemble  average 
of  (7)  in  terms  of  the  pair-distribution  function  pp(R„), 

<F,(f,  k')>,  =/(r,  k')  +  \pjd  Rs  p(R,J 

xj/(f',ft)</-s(fc,k'))s,e"k  (11) 

Proceeding  as  for  (8:13),  we  approximate  (Ft)u  by  (Fs)s, 
and  use  briefer  notation  F{ r,  k')  =  F^,  F(ic,k')  =  Fc0,  etc., 
to  obtain 

F«  =/o  +  \  pJrfR/KR)J/r^o«*-  - k|  R  (12) 

In  terms  of  the  Fourier  transform  of  the  pair  correlation 
function 

P  [K  ]  =  p  J /?(R)e,X  H£/  R,  K  =  k  (r  -  k)-(xx  +  yy) 

(13) 

and  simple  contours 

i  rn  rn/1  •  ioc 

—  d<Pc  ddc  sin  6C , 

llT  J  -  ir  JO 

we  obtain  the  form  (6:37)  developed  originally  for  striated 
surfaces 

P hi  —  f«i  +  8/„.  F ’o, 

S  =  1  lim  f  «  *  P  [Kc }  =  \-*P  [Kc  ]  +  S„ 

2  f  -o  jc  2 

(14) 

with .  ^  as  in  (4).  [See  (8: 1 3)— (8: 16)  for  discussion  of  the  con¬ 
vergence  factor  and  other  aspects.]  We  showed  that  Jc  oper¬ 
ates  on  a  function  of  the  real  variables  q>c  and  sin  6C ;  because 
we  assume  that  p(R)  has  the  same  inversion  and  reflection 
symmetries  as  an  ellipse,  P  [A"  ]  is  real,  and  it  follows  that  y  P 
and  S,  are  real  and  imaginary  operators,  respectively. 

Introducing  the  Fourier  transform  of  the  total  correla¬ 
tion  function  p(  R)  —  1, 

&[K]=pj[p(R)-l]e**dR,  (15) 

we  make  the  specular  contribution  explicit  to  obtain  form 
(6:41), 

FfO  =/o(l  +  C/v,l  + 

sm  =  yreP[Kc\  +  s,.  (16) 

If  the  effects  of  S101  are  negligible,  then  /"and  the  correspond¬ 
ing  impedance  reduce  to  the  sparse-gas  approximations' 
F-F*  and  Z~Z°,  with 

F%  =/o(l  +  CFSo)  =/o/(l  -  Cf), 

Z°  =  Cf[V,  k')  —  Cf,  (17) 

as  discussed  and  applied  earlier  in  detail.1-2  For  the  general 
case  (16),  we  reduce  F  by  using  the  results  of  the  analogous 
development  of  G  for  the  monolayer.8  As  before,  we  first 


introduce  the  energy  functions  to  provide  motivation,  but  in 
view  of  the  details  that  are  available,8  the  present  develop¬ 
ment  is  brief. 

II.  AVERAGE  ENERGY  FUNCTIONS 

We  write  the  average  energy  density  as 
<  |  |2>  =  |  ( >  |2  +  V,  with  |  <  V )  |2  obtained  from  (9)  as  the 

coherent  component.  The  variance  V  =  (\°2s  —  {fy)  | 2 ) ,  the 
incoherent  or  fluctuation  component,  follows  from  (6)  by 
proceeding  as  for  (8:19).  In  particular,  forz~  oo  we  write  V 
as  an  integral  over  all  back  space  directions  [\6  |  <w/2), 

V~  pc,J  7T[K  ]  \F(r,  k')|2  sec  6  =pj  Q*  sec  6- 


yP~[K  ]  =  \  +  &>[K)  =  1  +  pj [  /KR)  -  1  )eK*d  R,  (18) 

where  7F\K  ]  is  the  structure  factor,  and pQ ^  is  the  multiple 
scattered  differential  cross  section  per  unit  area. 

For  the  corresponding  average  energy  flux,  we  have 

J±  =  Re(V'*  W  ±  /ik)  =  Jc±  +1, 

Jc±  =k'±Re^'<Y'a±>*(k  +  k')+  |<Y'l>|2k, 

I~pJ  Qk)  r  sec  6.  (19) 

The  first  term  of  the  coherent  component  Jr  is  the  incident 
flux,  the  second  arises  from  interference  of  the  incident  and 
reflected  waves,  and  the  third  is  the  specularly  reflected  flux. 
The  net  returned  flux  is 

Jr=*k  +  I,  /J  =  |<V'r)|2  =  |1  +  2C/'|2=  ||±J-|2 

(20) 

corresponding  to  coherent  reflection  plus  incoherent  scat¬ 
tering.  The  components  of  (19)  normal  to  the  base  plane 
equal 

i-Jc  =  (1  +  R  )cos  60,  z-I  =  poa\-#  W\F^  |2  =p  y  „ 

(21) 

with  p  JP ,  as  the  average  incoherent  scattering  cross  section 
per  unit  area.  The  corresponding  average  absorption  cross 
section  per  unit  area  is  the  average  inward  flux 
—  z-J  =  p  y  a .  Conservation  of  energy  requires,  z-J  =  z-(  Jc 
-I- 1),  from  which' 

i  =  r  +Pysec0o,  y  =  ya  +  y„  (22) 

i.e.,  the  incident  flux  density  equals  the  flux  coherently  re¬ 
flected  plus  that  absorbed  and  incoherently  scattered  by  the 
area  of  distribution  irradiated  by  unit  area  of  the  incident 
wave.  From  the  development  following  (8:23),  we  may  repre¬ 
sent 

yB  =  -  J  Re('P*V'P/ik),-dA. 

as  an  integral  over  the  surface  of  the  boss,  or  as  one-half  of 
the  result  of  integrating  over  the  boss  plus  its  image. 

Substituting  R  =  |l+2CFoo|2  in  (22),  and  using 
4C  cos  QfJ  p  =  cr0,  we  obtain 
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-  Re  Foo  =  C  IFool2  +  (ya  +  y s)/a0 

=  c  |fjJ  +  ynr[K  \  |^|2  +  yjog, 

(23) 

which  for  lossless  scatterers  reduces  to  (6:65).  See  (8:24)ff  for 
the  detailed  procedure  for  the  more  general  problem  of  the 
nonsymmetrical  monolayer.  The  essential  feature  of  the  pro¬ 
cedure  for  present  purposes  is  that  the  energy  relation  (23) 
corresponds  to /'of  ( 1 6),  and  that  the  same  procedure  may  be 
applied  to  a  modified  problem  shorn  of  specular  losses. 

III.  SCATTERING  AMPLITUDE  REDUCTIONS 

The  procedure  (8:24)ff  for  obtaining  the  energy  rela¬ 
tions  (23)  applied  to  a  modified  scattering  problem  specified 
by 

A  =/o  +  S /«,,  s(0’  =  [K  \  +  s„  (24) 

i.e.,  by  the  form  (16)  without  the  specular  term,  yields  the 
corresponding  form  of  (23).  Thus  for  lossless  scatterers, 

-Re/00  =  \^7r[K}\/r0\2-  (25) 

In  terms  of  /we  reduce  (16)  and  the  corresponding  imped¬ 
ance  Z  =  CF /[\  4-  CF)  to 

FfQ  =/o(l  +  CF0 0)  =^/(l  -  C/), 

Z  =  C/=  C/oo  =  C/(k,  k'),  (26) 

i.e.,  the  same  forms  as  in  (17)  discussed  earlier,1  but  with  / 
replaced  by  the  modified  amplitude  / 

To  analyze  /  and  display  the  implicit  physics  we  ex¬ 
press  /in  terms  of  the  amplitude /'  for  a  radiationless  boss,4,3 

/,<,=/;<>  +HV|r/o.  (27) 

such  that  for  the  lossless  case,  f'Qy  + /|J  =  0.  The  leading 
term  of  (27)  includes  phase  and  absorption  effects  but  no 
radiative  losses.  For  a  lossless  isolated  boss, 

-Re/00  =  M'I/*>I2.  (28) 

with  the  special  case  for  Re/,,,  corresponding  to  (4)  for  oa 
=  0. 

Writing  i?[K]  =  %  [K  ]  —  1  in  S(t”  of  (24),  we  express/7 
in  terms  of /'  as 

/*,  =/*  + 1 K  r  +  %]/*/«■  (29) 

Suppressing  S,  by  means  of  the  radiationless  distribution 
amplitude, 

/*=/*+ 8. /Jr /5>.  (30) 

such  that  for  the  lossless  case,  Re /£,  =  0,  we  obtain 

Ao=fr0+^^lK]KA  (3D 

as  the  analog  of  the  isolated-boss  form  (27). 

From  Fm  =  //(l  -  C /)  =  ( /-  C  |/|2)/|  1  -  C/|2 
of  (26),  we  construct  —  Re  F00  =  ( —  Re/+ C\ /\2)/ 

1 1  —  C  / 12.  Using  (25),  we  obtain 

-ReF00  =  C|F00|2-(-l^^[^l|^|2  (32) 

corresponding  to  (23)  for  lossless  bosses.  In 
-  Re/op  —  pf'l/ol2  of  (28),  the  term  —  Re/*,  represents 
the  energy  loss  from  the  incident  and  image  wave  via  inter¬ 
ference  with  the  radiated  wave,  and  J^|/|2  shows  it  is  bal¬ 
anced  by  scattering.  In  (32),  the  interference  loss  is  balanced 
by  specular  reflection  and  fluctuation  scattering. 
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The  coherent  power  reflection  coefficient  equals 
R=  1  +z  2  -  i  4 Re Z  pfo  +g,) sec6»0 

1-Z  [1  —  Z  I2  |1-Z|2 

Z  =  C/oo,  a,  =  V\  /„  |2.  (33) 

The  corresponding  differentia]  scattering  cross  section  per 
unit  area  is  approximated  by 

P&o=P$,o/|l-Z|2,  qM  =  ariATJc.lAol2; 

c,  =  2/irk,  l/k  2,  (34) 

such  that  a,  =  fq with 


for  ar  isolated  boss.  The  phase  change  (0 )  introduced  by  the 
distribution  of  bosses  is  given  by 

tan  ©  =  2  Im  Z/(l  —  |Z  |2),  (35) 

where  0+  and  ir  +  <9  _  are  the  total  phase  changes  on  reflec¬ 
tion  of  ±  <f> '  from  distributions  on  rigid  (  4- )  and  free  (  —  ) 
base  planes.  These  forms  and  C  =  p/ky0  or pv/k  2y0  are  the 
same  as  before,1  and  the  dependence  of  / ±  on  y0  =  cos  0o 
and  y  =  cos  6  near  grazing  are  determined  by  the  corre¬ 
sponding  behavior  of f±  discussed  for  (1:11).  Thus /_(?,  k') 
vanishes  and  /+  does  not,  as  either  r  or  k'  approaches  graz¬ 
ing;  consequently,  Z+— *-oo  and  Z_— *0,  so  that  (IP*)/ 
(  ±  <t>  'H '  —  1.  and  Q—> 0. 

The  operation  over  S,  in  (30)  is  negligible  for  large  kb 
for  the  present  problem  (which  excludes  complex  60,  as  well 
as  multimode  periodic  cases),  and  also  for  small  kb  provided 
that  boss  width  (a)  is  small  compared  to  minimum  separation 
b.  (If  a/b  is  not  small,  S,  may  introduce  ^-independent  fac¬ 
tors  arising  from  multipole  coupling.8)  Neglecting  S,  corre¬ 
sponds  to 

Ao~f\o+^r-f\r/r2,  (36) 

which  differs  from  the  isolated  scattering  amplitude  /,0  in 
that  the  radiation  integral  is  modulated  by  the  structure  fac¬ 
tor  —  1  +  & .  As  discussed  before,8  the  term  1  corre¬ 
sponds  to  radiation  in  isolation,  to  correlations  among 
pairs,  and  7^  represents  fluctuation  scattering;  see  (8:64) 
and  (8:77). 

For  small  Im  f> Re /,  or  small  i^'and  small  absorption, 
we  iterate  (37)  and  neglect  absorption  in  the  quadratic  terms 
to  obtain 

A,  -foo  +  =/»  - 

Im/oo  -  K  +  os)/<70. 

(37) 

To  this  approximation,  /’differs  from  /only  by  the  presence 
of  the  structure  factor  7F\K  ]  in  the  radiation  integral.  Writ¬ 
ing 

Z  =  C  /go  =  «Z]  +  Zj, 

Z,  sC  Im/oo  =pcr0  lmf/4 y0, 

Z2=r -pier.  +  a,)/4y0,  (38) 

we  simplify  (33H35)  for  various  practical  purposes  by 
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R  j  _  P\°*  +  a,  )sec  d0 

~  |l-Zp  ’ 

„  p7T[K\q^  _  0  „ 

fQ=  irTz,-  ■  '“Y=z"  (M| 

where  we  may  approximate  1 1  —  Z  | 2  by  1  —  2Z2  +  Z] 
si  +  Z 2 .  The  factor  sec2  6  in  Z2  insures  that  /?  +— >1  as 
00—-ir/2,  and  also  relates  to  the  pseudo-Brewster  minimum 
in  R  + .  The  present  Q  may  be  used  to  isolate  77\K  ]  versus  r 
from  measurements  essentially  as  in  x-ray  diffraction  by  li¬ 
quids,  and  the  resulting  Tf  —  1  =  &  may  then  be  Fourier 
transformed  numerically  to  determine  p(R)  —  1.  For  the  di¬ 
rect  problem  we  use  W'\K  ]  as  in  (8:76)  and  (8:77),  or  numeri¬ 
cal  results  obtained  from  other  approximations10  of  p. 

For  cylindrical  bosses,  except  for  to  =  pbx  1,  it  follows 
from  the  Zemike-Prins  p(x),  that6-® 

=  1  +  2w  f  (p  —  l)cos  u  du  =  (1  +  3)~\ 

Jo 

3  _  to  sin  w2(  1  —  cos  JT) 

~T~{1  -to)JT  +  (1  —  u»)2^2  ’ 

JV  =  kb  (sin  6  —  sin  ff0).  (40) 

If  kb  is  small,  or  if  G0,  then 

77  \V)  =.  77  +  BJT 2  +  tfpr4),  7T  =  (1  -  to)2, 

B  =  77  w\(\  -  \w),  (41) 

so  that  77  increases  quadratically  in  3f  from  its  mini¬ 
mum  77  =  77  (0).  The  first  local  maximum  of  77\,7(/')  of  (40) 
with  increasing  is  the  largest,  and  for  large  JFf  1  —  to), 

»m~ i-,i»  ,42, 

1  -  w  jr 

approaches  unity  with  oscillations  decreasing  as  3f~ '. 

For  aligned  bounded  bosses  and  a  circular  exclusion 
curve  of  radius  b,  corresponding  to  hard  statistical  mechan¬ 
ics  disks  of  diameter  b,  we  have 

7/  '(JK)  =  1  +  8to  f  {p  —  l)70(JF'n)u  du, 

Jo 

w=pirb2/ 4,  (43) 

/kb  I2  =  (a  -  a„):  +  ( /?  -  0of 

=  sin2  6  +  sin2  0„  —  2  sin  6  sin<?0  cos(g?  —  <p0), 

where  Jf2  =  b  2K-(ix  +  yy)-K  with  K  =  k  (f  —  k).  In  gen¬ 
eral,  computations  for  7P~(Jr/')  can  be  based  on  numerical 
results  for  p  obtained  from  integral  equation  approxima¬ 
tions.  10  If  kb  is  small,  or  if  r~  k,  we  may  use  the  virial  expan¬ 
sion  ofp  to  develop  W(W)  in  powers  of  to  and  JT2, 

5r-(jr)=:  1  -4u;  +  6.6159to2  +  {JT2w{  1  -  5.9746 to) 
s  7T  +  3T2B,  (44) 

where  the  final  form  is  the  analog  of  (41).  From  the  scaled 
particle12  approximate  equation  of  state,  we  showed 

7T  =  5r(0)~(l  -  io)3/(l  +  to)  =  1  -  4ur  +  Iw2  + 

For  large  ( 1  —  w)JT,  we  integrate  (43)  by  parts,  and  use 
y,(jc)~(2/jrx),/2  sin(jt  —  rr/4)  and  the  scaled  particle  approx¬ 
imation  for  p(b )  to  obtain 


»  ir)~\  - 


\  _  4w(2  —  w)(2/it)U2  sin (JT  -  ir/ 4) 

,  (1  —  w)2  JT3'2  '  '  ’ 

which  approaches  unity  more  rapidly  with  increasing  JT 
than  (42).  For  elliptically  symmetric  statistics  specified  by  an 
exclusion  ellipse  with  principal  semidiameters  b2  and  b}  (cor¬ 
responding  to  aligned  hard  elliptic  disks  with  principal  di¬ 
ameters  b2  and  f>3),  from  geometrical  consideration,®  we 
work  with  (43)— (45)  in  terms  of 

3F2/k  2  =  b\(a  —  a0)2  +  b 2 ( (3  -  0O)2,  to  =  ir pb2bi/ 4, 

(46) 

where  Jf2  =  K-(b  2xx  +  b  2yy)-K.  We  apply  existing  results 
for  circular  disks p(R  )  =p{  u  j  given  in  terms  of  u  as  radial 
distance  divided  by  disk  diameter,  by  regarding  u  as  the  size 
parameter  for  a  set  of  cocentered  similar  ellipses;  u  =  1 
specifies  the  exclusion  ellipse. 

For  small-spaced  bosses  (small  kb ),  and  small  a/b,  to 
lowest  orders  in  k  for  the  numerators  and  denominators  in 
(39),  we  have 

B~1  _  p[a°  + 

1+Z2  ’ 

,47, 

1  —  Z  f  4  cos  dQ 

where,  except  for  pressure  release  bosses  on  a  rigid  base,  Z,  is 
of  order  k,  and  er,  and  q  are  of  order  k  3  for  cylinders  and  k 4 
for  bounded  bosses.  To  this  approximation,  the  greater  gen¬ 
erality  of  the  present  distribution  compared  with  the  sparse 
gas  case,1  appears  simply  as  a  multiplicative  packing  factor: 
thus  the  earlier  results  can  be  carried  over  by  inspection. 

The  packing  factor 


»  =l+pj  MR)-1] 


decreases  the  scattering  effects  as  the  packing  density  w  in¬ 
creases.  If  there  are  n  bosses  in  a  central  region  A,  such  that 
p  =  {n)/A,  then 

=  [<n2>  -  {n)2\/ (n)  =  <1  n  -  v n))2)/(n >  (49) 

is  the  variance  (fluctuation)  in  the  number  of  bosses  in  A.  If 
there  are  no  fluctuations,  then  Q  and  y,  =  7/  a,  /( 1  +  Z  ] ) 
vanish,  and  R  <  1  corresponds  solely  to  absorption.  For  hard 
statistical  mechanics  particles  governed  by  an  equation  of 
state  E,  we  used  the  theorem9  ( dE/op)~ 1  =  7/  and  the 
scaled  particle12  approximations  frv  L  toobtain"  ?/  ’ as  sim¬ 
ple  rational  functions  of  the  par'.ing  density  w: 

ar0  =  (i  -u>),  7?',  = . l  -  u>)2, 

ar2  =  (l-u;)7(l  +  u;).  (50) 

Here  ,  is  the  function  in  (41),  77\ is  mentioned  after  (44), 
and  W0  is  the  analog  for  a  random  lattice  gas  corresponding 
to  uncorrelated  space-occupying  particles.  The  upper  bound 
on  7E~ o  and  7/  ,  is  unity;  for  we  used"  to  5  0.84,  an 
experimental  value  for  the  densest  random  packing  of  identi¬ 
cal  circular  disks  on  a  plane,  which  approximates  the  mean 
of  the  square  (0.785)  and  hexagonal  (0.907)  close  packed  val¬ 
ues. 

The  fluctuation  effects  are  determined  by  ( n )  W~,  or  by 
the  normalized  variance,  S  =  w7E\  From  (50),  we  see  that 
S  (to)  is  small  for  to  near  its  lower  and  upper  bounds,  and  that 


89 


J.  Acouat  Soc.  Am.,  Vol.  73,  No.  1,  January  1983 


Victor  Twersky:  Scattering  by  rough  surfaces 


89 


5  has  a  maximum  =  S  (wn )  corresponding  to  w  =  wn  for 
which  dwS  =  S'  ~  0.  As  discussed  before  in  a  different" 
context, 

V,  =  4>  a'(y ,  = 

S,n  =  AsO.148,  u;ln  =  ); 

•S^sO.0856,  Wi, -0.215.  (51) 

To  consider  extrema  of  R  and pQ  as  functions  of  w,  we  write 
(47)  as 


R  =  1  - 


wF + sb 
1  +w2A2 


1  f  u/A  2  ’ 


and  assume  F<5. 

For  moderate  #0  we  drop  A  2.  Then,  if  there  is  no  absorp¬ 
tion,  R  —  1  —  SB  has  a  minimum  and  pQ  —  SB '  has  a  maxi¬ 
mum. 


R^\-SnB,  pQr~SrB\  (53) 

at  the  values  w,,  of  (51)  corresponding  to  maximum  fluctu¬ 
ations.  For  non  vanishing,  small  r  /B,  the  extrema  are  shifted 
to  wA  >  w  .  To  first  order  in  F,  in  terms  of  5 "  =  d 2  S, 

R^l-(Wnr+SnB),  w^~wn+F/B\S''\,  (54) 

where  1/|5”|  =  J,  1,  0.453,  for  the  sequence  in  (51).  For 
small  w,  scattering  losses  dominate  in  R,  but  with  increasing 
w  absorption  dominates.  On  the  other  hand  if  F  /B  is  not 
small  then  there  is  in  general  no  minimum  of  R  with  vari¬ 
ation  of  w,  and  r  dominates  the  loss  effects  for  all  w. 

Similarly  if  we  retain  A  2  in  (52)  for  F  =  0,  we  obtain 
frequency  dependent  values  wa  for  the  extrema  correspond¬ 
ing  to  vanishing  Rl~R'-,  +{ww  —  w,, )R  To  lowest  or¬ 
ders. 


R„  s  1  -  S  B  ( 1  +  u&A  2),  pQ„  ^SnB '/( 1  +  w2A  2), 
wIO  -  w,  ~-R  'JR  -  Sr,w,2A  2/|S"|.  (55) 

Thus  the  extrema  occur  at  smaller  packings  ww  <wn,  and 
the  shift  is  proportional  to  A  2.  Extrema  with  both  T  and  A  2 
retained  correspond  to  w  such  that  in  terms  of  R  (u»n )  =  R  : 
and  Q(wn)  =  Q,  and  of  (52),  we  have 

Rj  ~R~t  Qj  ~  Q  , 

wj  -  2 wnSnA  2/|5"|  +  F/B  |S"|  (56) 

when  wA  may  be  greater  or  less  than  wn. 


IV.  SCATTERING  COEFFICIENT  REDUCTIONS 


In  Ref.  8,  we  expanded  the  scattering  amplitude  g  of  an 
arbitrary  obstacle  and  the  multiple  scattering  amplitude  G 
for  a  general  monolayer  as  multipole  series  in  terms  of  corre¬ 
sponding  scattering  coefficients  a  and  A .  From  G  [g]  of  (8: 16), 
we  obtained  an  algebraic  system  A  [a]  expressing  A  in  terms 
of  a  and  the  distribution  integrals  (continuum  analogs  of 
lattice  sums  determined  by  p ).  The  same  procedure  applied 
toF  [/]  of  ( 1 4),  or  more  directly  to  Fof  (5)  in  terms  of  G  for  the 
symmetrical  monolayer,  involves  the  same  A ’s  and  ^s. 


'lii us,  from  (5),  for  semicylinders 

F%  =  Gm  ±  G,n  =  £  An  [e‘"e  ±  <?'"<"  -  *•],  (57) 

n  »  0 

with  A„  ( 60 )  as  in  (8:57)ff  in  terms  of  a„  and  „ ,  see  (8:59)— 
(8:70)  for  in  terms  of p(x),  or  in  terms  of  W\K  ].  Similarly 

for  bounded  bosses, 

F%=  £  £  A:[Y:(r)+Y”[r')) 

<1  =  0/1=  —  m 

=  ^A:Y”(i)[  i  +  (-ir— j,  (58) 

nm 

with  A  "(r0)  as  in  (8:72)ff  in  terms  of  a”  and  see  (8:74)- 

(8:89)  for  38J.  All  of  the  earlier  results  for  the  symmetrical 
monolayer8  can  be  applied  to  the  present  problem.  We  illus¬ 
trate  the  procedure  by  considering  semi-ellipses  and  hemi- 
ellipsoids  with  principal  diameters  d,  small  compared  to  A. 
In  the  course  of  the  development,  we  provide  corresponding 
forms  of and  /"to  be  used  in  (26)ff. 

For  ellipsoids,  or  for  other  obstacles  having  the  same 
reflection  and  inversion  symmetries,  we  retain  only  the  mon¬ 
opole  and  dipole  terms  and  write  the  scattering  amplitude  as 

8, o  =  a0  +  £  a.Y.fr, 

i  =  1 

a0  =  ai/(l -a'0),  a,  -  a',/{\  -  a'/n),  (59) 

where  a'0  and  a’  are  the  corresponding  coefficients  of  g'  for 
the  radiationless  problem.  Here  n  =  2  or  3  for  cylinders  or 
bounded  obstacles,  and  yi  and  are  the  direction  cosines  of 
r  and  k  =  r0  with  respect  to  z,  x,  y  [see  (8:90)  for  arbitrary 
alignment).  For  a  boss,  from  (3), 

/*  =2a,rltf, 

/*  =  2(a0  +  a2y2y°2  +  a, ntf )-  (60) 

with  similar  forms  in  terms  of  a'  for  the  corresponding  /'. 
For  the  distribution  of  bosses, 

F  k>  =  ^AiYit 

F  =  2(A0  +  A2y2  +  A3y3),  (61) 

with 

<4 1  =  ailY?  +  38  uA,),  A,  =  aly?  4-  3F ,jAt\ 

V  i  / 

jr„  =  . ^,=2S)^,  y0=  1,  (62) 

where  randy  equals  0,  2,  3  in  the  form  for  A,,  and  the  38, j  are 
given  in  terms  ofJFn  and  38^  in  (8:93)  and  (8:94). 

For  A  |,  we  have  directly 

F =  2.c/ly,y?  =  2.^,  cos  ^cos  0O.  (63) 

The  three  coupled  equations  for  A(),  A2,  and  A  ,  are  reduced 
after  (8:95).  In  particular,  if  <p0  =  0,  so  that  =  0  and 
=  sin  60 ,  then 


A o  —  ^oll  +  '^02^2)’  A2  —  J 2(y?  4-  JF 2o4o)i  sf  t  /( 1  —  at38  „ ), 

p  +  _  j  ‘^0  yjy 2  ~F  )  _  2  ao  028282  —  00^2 [^22  +  ^00 7282  ~  ^ zniX;  4-  y^)] 

1  -  ^0^2(^02)2  ~  1  -  -  a2jr22  +  flo«2[^ oo^22  -  (^o:)2] 
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This  result  holds  for  all  k  for  spheroids  with  symmetry  axis 
along  z,  and  for  the  two-dimensional  problem  (y2  =  sin  6 )  of 
elliptic  cylinders. 

More  generally,  we  reduce  the  coefficients  to  determine 
the  modified  amplitude  /  to  use  in  (33)ff.  To  isolate  the 
specular  contributions,  we  decompose  the  distribution  inte¬ 
grals  as  in  (8:105), 

=  2 C-fifj  -  m,S,t  +  JT'j, 

jr'u  =^;,+i'v  <65> 

where  mQ  —  1,  and  m,  =  1/n  with  n  =  2  or  3  if  i^O.  Substi¬ 
tuting  into  the  initial  form  of  A,  in  (62),  and  using  F w 
=  2A,rf  of  (6 1 ),  we  obtain 

+  [2C(tf):+  -  1/«]^,| 

=  «'|7'?|1  +  00  )> 

«,■=«»,/[  1 —0,(^1,  —  1/n)].  (66) 

Thus,  from  FM  =  /H) ( 1  4  CF(10)  of  (26),  and  F^  ~  2A  ,y,  of 
(61),  we  have 

=  2«iy,r?. 

«,  =  a|/(l  -  a\  Ft )  =  aj'/(l  -  a"  fn),  (67) 

whereaj  —at/(  1  +  £7,/n)andn;'  =  a\/{\  —  ia /',,)arethe 
corresponding  coefficients  of  /'  and /". 

For  F  * ,  we  introduce 

a  =  a,xx  4-  a .yy,  A  =  a,x  4  a, y, 

*  = 

¥  =  /f  ,  ,xx  +  ¥ , ,|xy  4  yx)  4-  2^  ,  ,yy,  (68) 

or  the  corresponding  matrices,  and  rewrite  the  system  in  (62) 
as 

Av  =  n0(l  +  &  (Hy4(,  4-  ¥-A), 

A  =  a-(k  4-  ¥An  4-  .¥-A).  (69) 

Using  (65),  and  compacting  F ^  —  F  as  in  (61),  we  write 
y,x  4-  y} y  as  k'  (or  k)  to  obtain 

Aq  =  ■'/  o  ( 1  4  CF+JT-A), 

A  =  ,y-[k'(l  4-  CF I  4-  ¥A„\. 

A  a  =  ou/(  1  —  o0.¥  no )  =  a,,  /( 1  —  o0  f  oo ), 

A  =  (!  -  i'-J'V  =  (/-  a "■/')-  ‘-a*.  (70) 

Here  a'Q  =  n,/(l  4-  a0)  and  a'  =  (I  4-  i /n)~  ‘-a  =  a'2xx  4-  a\ 
yy  are  the  coefficients  of/',  and  a",  =a'0/{\  —  ia'0.  t  (K))  and 
a"  =  (/  —  /i\  are  the  coefficients  of/’.  Eliminating 

the  cross  terms, 

A o  =  f«o  +  )U  +  CF), 

A  =  («  4  «-k  )( 1  4-  CF ), 

»0  =  d-  ,yjjr..5r.jn  =  .t/-(i  4- 
=  jr.,rfat,  =  .tf-JTa,,.  (71) 

From  (26)  and  (71), 

/kS  =  2[<*0  4-  «•(?  4-  k  )  4-  r-a-k  J,  (72) 

where  only  the  £  and  y  components  of  ?  and  k'  arise.  We  may 
also  obtain  (72)  by  using  /,  =  2(a0  4-  Fa/,)  and  /*, 
=  2(Uo  4-  ?•»)  in  (24)  to  obtain  u  =  a0  4  a- k'  and  u  —  a 
4-  «-k'.  See  analogous  discussion  of  y  following  (8: 109).  Spe¬ 


cializing  to  <p0  =  0,  reduces  (72)  to  the  form  (64)  in  terms  of 
o ,  ■¥ 2  =  /(I  —  F  22 ),  and  ■¥  t 02 . 

We  may  apply  the  above  to  Rayleigh’s  results15  for  el¬ 
liptic  cylinders  and  ellipsoids  of  area  or  volume  F",  specified 
by  two  relative  parameters  C  and  B  (compressibility  and  in¬ 
verse  mass  density  in  the  simplest  cases,  but  complex  in  gen¬ 
eral).  We  have 

a'0  ~iFF,  f€  =  C-  1; 

a;  ~  _  ,7W(  14 -  j&q,)  =  -  iF 93,',  .3  =  B  -  1, 

r  =  k2?"7 4,  k2F7 4jt;  Im  C>0,  lmi?<0.  (73) 

These  results  are  correct  to  /  (A "),  and  the  next  terms  are 
¥{k  "  '  2).  Here  q,  is  the  standard  depolarization  integral1'’ 
written  in  terms  of  the  principal  diameters  d,. 


q,  =  d,d,d,  f  j  2(d 2  4-  x) [{d  {  4-  x)(d 2  +  *) 

JO 

Xld;  +x)]',2j  1  dx,  j?  q,  =  1. 


For  cylinders,  d,— *oc  and  =  0, 


1  -?>  = 


1 


d ,  4-  d,  1  4  r 


(74) 


(75) 


with  q |  =  <jr2  =  (  for  the  circle.  For  spheroids  with  d,  as  the 
symmetry  axis,  <?,  4  2q:  =  1.  If  r  =  dx/d,  <  1,  the  oblate 
case,  then 


q{  = 


r cos  1  r  ^ 
"(I  _  r=)*~/ 


3 


<<h  <  •. 


0<^2  < 


3’ 


(76) 


such  that  for  near  disks,  q,  x  1  —  nr/2  4  2 r.  If  r>  1,  the 
prolate  case. 


1  j  r  cosh  1  r  , 
q'~  ^_i  l^-l)'  2 

0<9.  <y.  y<?:<y.  (77) 


such  that  for  near  needles,  9,~(ln  2r—  1  )/r.  If  d2  is  the 
symmetry  axis,  then  q2  +  2q,  =  \,  and  we  interchange  d , 
and  d2. 

If  r  =  1 ,  then  q,  =\  for  the  sphere.  For  near  spheres, 
r~  1,  and  (76)  and  (77)  give  9,  —  4 [r  —  1  )/l  5.  More  gener¬ 
ally,  for  a  near  sphere,  d,  =  d(l  4  <5,)  with 5,  ~0.  From (74) to 
first  order  in  5,  we  have  q,  =  J  4  (  —  4<5,  4  25,  4  25,);  simi¬ 
larly  for  q2  and  q„  so  that  2q,  =  1.  For  a  volume  preserving 
deformation  (d,d,d,  =  d  ’),  we  require  25,  =  0,  and  conse¬ 
quently  q,  =  (|  —  25, /5.  These  points  are  stressed  to  facili¬ 
tate  geometrical  interpretation  of  the  multipole  coupling  ef¬ 
fects  we  consider  subsequently. 

The  form  in  (73)  for  Oo  holds  for  arbitrary  shapes;  if 
then  a,'  ~  —  iF./)  is  similarly  shape  independent.  For 
a  rigid  obstacle,  we  let  C  and  B  equal  zero  in  (73)  to  obtain 


a!,  s:  -  ir,  a’zziT /(I  -  q,\.  (78) 

For  pressure  release  obstacles  (e.g.,  gas  bubbles  for  under¬ 
water  sound) 
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flo2  ~  -  iir/ 24  floi  ~  -  t'k/  a,'  =  ~ir/qit 
f0  =  ln[8/kc(d,  +  d2)],  c=  1.781; 

'=  (J[”  [\d\  +  x)(d\  +  x)(d\  +x)]-l/2dxy\  (79) 

The  dipole  terms  a(,  obtained  by  letting  B— *00  in  (73),  are 
negligible  in  g  but  dominate  in /“.  The  values  of  /(the  elec¬ 
trostatic  capacity)  for  the  spheroid  equal 

/2/d2  —  (1  —  r*),/2/cos_’  r  for  r<  1, 

/2/d ,  =  (r2  -  l)l/2/cosh~" 1  r  for  r>  1.  (80) 


These  reduce  to  /2/d  =  1  for  the  sphere,  to  /2/d2~(2/ 
rr)(  1  4-  2/VV)  for  near  disks,  and  to  /2/d,  s  1/ln  2/-  for  near 
needles. 

More  generally  than  in  (73),  we  use  a  dyadic  parameter 
B  and  incorporate  the  shape  integrals  q,  in  a  dyadic  q.  For 
present  purposes,  we  assume  that  the  principal  axis  of  all 
dyadics  are  aligned  with  the  Cartesian  axes  of  the  exclusion 
ellipse  and  its  normal,  and  replace  28  by  =  Bi  —  1.  See 
Ref.  8  for  the  general  case. 

For  small  spaced  scatterers,  we  use 

fr-TITmA,,  /ysvM,.  («U 

where  7f'  is  the  k-independent  packing  factor.  The  -  ru  for 
k:= 0,  reduce  to.  1  oo—^o.  7i°/kb2l  and for 
i#0;  except  for  310  which  depends  on  In  k,  the  3l’s  are  fc 
independent. 

For  such  cases  (for  /  =  0,  1,  2,  3), 

<  „  a'i  a'i 

a,  ss  — - - ,  a  2: - =  — —  .  (82) 

1  —  a"  78  m,  \  —  ia’nu  2Zt 

Except  for  the  monopole  terms  for  pressure  release  scat¬ 
terers  (which  we  consider  separately)  all  coefficients  are  of 
order  k".  The  dominant  terms  in  (81)  for  cylinders,  from 
(8:70)  and  (8:93)  are  given  by  (kb ) 


in  terms  of  the  Zemike-Prins  p.  If  w  is  near  zero,  then 
/2s2u)(1  4-  w  0.307);  if  w  is  near  unity,  /2=;f(2)  =  it2/ 
6-1.645. 

For  bounded  obstacles  and  elliptically  symmetric  sta¬ 
tistics  determined  by  /  =  b3/b2  (the  ratio  of  the  principal  di¬ 
ameters  of  the  exclusion  ellipse),  from  (8:87)  and  (8:94)  in 
terms  of  (kb2)2  k~u  — ►9?,, , 

H„=j*S  =  I^^rP-du, 

J 1  u 

[*”}  =  -  i*;  ±  m = w>2  ±  wi 

=  J/2(M°±3A/2), 

M°2=(2/irt)E(v), 

M\  =  (2/jr3fv)[(2  -  v)E(v)  -  2(1  -  v)ff(v)], 

v  =  1  -  r 2,  (84) 


where 


E(y\  = 


(1  —  vsin2  r)xn  dr 


and 

rv/  2 

/f(v)  =  J  (1  —  v  sin2  r)~'n  dr 

are  the  complete  elliptic  integrals,  and  p  is  the  two-dimen¬ 
sional  radial  distribution  function.10  If  w  is  small, 
/2s8u;(l  4-0.4157u;).  See  discussion  of  (8:89)  for  illustra¬ 
tions  for  t  =  2  and  4,  and  for  comparisons  with  analogs  for 
the  rectangular  periodic  lattice.5 

For  all  such  cases,  3> ^  =  1  +  0(k  2)z;  1,  and  the  mono¬ 
poles  satisfy 

<io  ~ao  —iPiS  ■  (85) 

The  dipoles  a"  =  a\/ 22> ,  differ  from  a'  by  the  k-independent 
factors 

3,  =  1  -  »  1  -  rs;7t„/(kb2r  =  1  + 

(86) 

For  cylindrical  dipoles,  from  (73)  and  (83), 

e,  =  e,  e2=  -e;  e=n2/2nb\  (87) 

where  =  ird3d2/4.  For  bounded  dipoles,  from  (73)  and 
(84) 

e ,  =  e,  e2=  —  J(e  +  /),  e3  =  -  \(e  -  /); 
e=n2M°2/Mrb\,  e’  =  in^l/lirbl,  (88) 

where  V  =  ird^dAi/b-  Both  sets  satisfy  =  0,  so  that  if 
we  introduce  new  depolarization  factors 

e,=?,+e„  x e,  =  1.  X/=°.  (89) 

we  reduce  o"  to  the  same  form  as  a', 


a;  -  ins;  -  trst, 
1+39,'e,  1  +  28, (q,  4-  e.) 


ires, 

i  +  a.Q, 


-ins,. 


(90) 


The  relations  (85)  and  (90)  specify  an  ellipsoid  having  the 
same  volume  and  physical  parameters  as  the  original  but  a 
different  shape.  From  the  discussion  following  (77),  the 
packing  effects  incorporated  in  Q ,  determine  the  shape.  The 
equivalent  ellipsoid  is  flatter  along  the  array  normal  (i  =  1 ) 
and  broader  in  the  base  plane,  the  elongation  being  greater 
along  the  smaller  diameter  of  the  exclusion  ellipse  (along  b2 
for  b2  <  £3). 

For  elliptic  cylinders  with  shapes  specified  by  r  =  d,/ 
d2  =  q2/qx,  the  equivalent  cylinders  correspond  to  r, 
=  d\/d’2  —  Q2/Qx  such  that 


rt  d\d2  \d‘,Y 

r  dr2d,  d, 


_  [  dl  l2  -  1  ~  f/?2 

ld2 .  \  +  e/qy 

\d1±d£r 
8  b2  2‘ 


(91) 


We  substitute  (85)  and  (90)  in  «,  of  (82),  and  write 
^  +  '‘Z (2i  and  ®,-  =  —  t®fl2l  in  terms  of  real 

components  with  #  2  and  58 2  positive  and  small  compared  to 
,  and  .  Keeping  only  the  leading  loss  terms,  we  have 


92  J.  Acoust.  Soc.  Am.,  Vd.  73,  No.  1,  January  1963 


Victor  T wersky:  Scattering  by  rough  surfaces  92 


(98) 


«o  =  «o/(  1  -asnns:irem-r[tigm  +  rar^f„], 

«.  =</d  -«ry/»)«  -£r»,„  -r(»m  +r^„) 


-jT33 


<nt 


r  /  rir@2  \ 

&?('- +-A4 


1  +  £$<,1,61  (i  +  ^ 

(92) 

Using  these  in  (67)  and  (68)  in  the  forms 

A  =  2(*0  +  ),  A  =  2*  inr1? .  (93) 

we  substitute  into  F  =  //( 1  —  Z )  with  Z  =  C/. 

Writing  y“  =  y,  for  brevity,  the  imaginary  part  of  Z  is 
given  by 

ErnZ±  =4^3-,?,; 

*Y i 

Sid  =  #111  -  s82<nrl  -  ®jmA  3i  =  —  $imA  (94) 

The  corresponding  real  part  involves  absorption  and  scatter¬ 
ing  terms 

-p{aa  +as) 


Re  Z  ^  =■ 


pk7~ 

2  r, 


rv±  -  P-2L2. 

8  y 


(95) 


O'lal  —  A)  +  '2?2|2|  y\  +  3?3(2(  7^ »  ril<J1  —  '®II2I?/1' 

.  .  .^2  ,  .  ^mH  ^ 

tfui  —  r5  in  H - 1 - .  rim  =— - 


In  a , ,  for  semicylinders  we  use  A: 3  and  n  —  2,  and  for  bound¬ 
ed  bosses,  k  A/ir  and  n  =  3. 

For  rigid  bosses  on  rigid  or  free  base  planes,  by  letting 
i>—*  —  1  and  ►  —  1  in  (92)— (95), 

>  +  _  ,  ,  ri  ri  . .  ri 

rim  -  —  1  +7; — ~  +  ~ — ~  >  rim  = 


riw  —  1  + 


(i  -Qi)  (i-e,) 

ri  ,  ri 


i  -  <?. 


(96) 


(i-&)«  u-&)2» 


t?l») 


Y\ 


(1  ~Q,fn 


The  corresponding/*  and  Z+  specify  reflection  for  the  ho¬ 
mogeneous  case  of  a  fully  rigid  rough  surface;/,  is  of  inter¬ 
est  for  underwater  sound  incident  on  the  submerged  part  of  a 
rigid  obstacle  on  the  sea  surface,  and  Z_  specifies  reflection 
for  the  corresponding  distribution. 

The  present  development  also  suffices  for  the  homogen¬ 
eous  pressure  release  surface,  free  bosses  on  a  free  base, 


5m  =  - r\/Qu  S(;  =  r\/Q\n. 


(97) 


For  various  practical  purposes,  we  would  use  3  +  of  (96)  and 
3  "  of  (97)  for  the  reflection  of  sound  from  the  sea  surface, 
3  +  for  airborne  sound,  and  3  ~  for  underwater  sound. 

The  one  remaining  case,  that  of  free  bosses  on  a  rigid 
base  plane  is  exceptional  in  that  the  monopole  as  in  (79)  is  not 
of  order  Ac".  Here/_=a„  is  appropriate  for  underwater 
sound  incident  on  a  gaseous  boss  adhering  to  a  rigid  bottom. 
The  effects  of  dense  packing  on  the  imaginary  coefficient  a'0 
as  determined  by  ~  Ao»  may  be  appreciable.  We  have 


ao  ~  ' 


=  4. 

kb2 


1  a0LV oo 

For  cylinders  yi^oo— 320, 

3?o=--/°,  1°  =  /0  In  - /,, 

ir  ckbu 

1°  =  2u;  H  (p  -  Din  du=Ii-  Is, 

Jo  ckbu 

I'o  =  2w  In  f  (p  -  1  )du  =  -  2w(  1  -  — ]ln  — , 
ckb  Jo  \  2  /  ckb 


(p  —  l)ln  u  du. 


(99) 


If  w  is  small,  —  2u>[l  +  w{2  In  2  —  1.25)];  if  w~  1,  then 
/•~In(2tr),  and  7° - ln(4rr/ckb ).  For  bounded  obstacles, 


*8  =  -  IoM 


0,  4  =  84 

Jo 


\]du, 


MZ=(2/nt]K(v),  v=l  -t~2.  (100) 

If  w  is  small,  then  70=:  —  8ui(l  —  0.64 13w).  For  the  rectan¬ 
gular  lattice,  we  replace  IqMZ  by  L  °  =:  —  2  ln(4rr/ 
cr )~  -  3.906  +  2  In  t. 

We  have 

^02=:1  +/°/4=  1  -  |/°|/4 

^03c=  1  +  (S/btfoMZ  =  1  -  (f/b2)\I0\M°0,  (101) 

where  only  is  k  independent.  The  corresponding  imped¬ 
ance  is 

Z*  =  CA  =  C2«0 

=  C2ao/{\  -  a;; 7T)  =  C2a'0/\£#  -  a’0  JT) 

=  C2a'0/{\  -iaW  -a’0W^,  (102) 

with  Op  as  in  (79).  For  cylindrical  bosses,  a'a  =  —  rir/24, 
Z+  =  -;2pAy,[2(4  +  I°)/ir  +  i»"],  (103) 

where  for  w~  1,4  +  1° ~\n[2b /ir\di  +  d2)]  becomes  inde¬ 
pendent  of  k,  and  W  =  (l  —  w)2~  0.  For  bounded  bosses,  a'n 
=  -ikf, 

z  =  -  (hr pf(  1  \ 


kyt  Vl  +  /IJ»i%/bi  +  tknrj 

—  ~  72tt pb2(  1 


(104) 


kyt  \b2/f+IoMZ  +  ikb2»J’ 

so  that  for  both  (103)  and  (104),  Im  Z  is  essentially  of  order  1/ 

k. 

We  generalize  the  monopole  development  to  include 
the  first  resonance  for  large  finite  #  =  C  —  1  ~  C,  by  using 


a6=  ■ 


ir<g 


24  1 

v  ’  k/ 


(105) 


(1  -F€¥\ 

The  denominator  is  of  the  form  1  +  ^  tf(k 2);  if  iS  is  moder¬ 
ate  in  size,  we  neglect  the  &\k 2)  term  to  obtain  (73),  and  if 
^—►oo  we  obtain  (79).  The  resonant  frequencies,  given  by 
TV  y  =  1,  equal 

k\„  =lir/7r<€  4,  4  vf/W  (106) 

and  the  corresponding  isolated  monopole  coefficient 
ao  —  ao/(l  —  fli)  reduce  to  a0A  =  -  1. 


J.  Acoust.  Soc.  Am.,  Vol.  73,  No.  1,  January  1983 


Victor  Twersky:  Scattering  by  rough  surfaces 


93 


93 


For  the  scatterer  in  the  array,  because  of  multipole  cou¬ 
pling,  the  resonances  occur  at  different  frequencies,  and  the 
resonant  magnitudes  are  higher.  From 

«0~flo/(l  -  a'0L  i  oo  -  a’o  ?/  ]  =  aSAl  ~  ao  (107) 

and 


a'0  _  ir€  _  rTV  _ 

1  —  OgL  I  oo  ~  1  —  - .  I'oo)  ~  \  —  k2/k\ 

(108) 


we  have 


jt2 

k\ 


1  -x 


_  rr%,  _  |/°| 


2/„ 


fo 


(109) 


andrtoA  =  -  1/^  ',  |«oa  I  >  Ka  I  =  1  •  The  reflection  coef¬ 
ficient /{=  |(1  +Z)/(  1  -Z)|2,  Z=C/~2C«0,  mk  =  kA 
reduces  to 


R,=l - — - -  =  1  -Pys  sec0o-  (HO) 

A  \»'+2Cl:  H 

With  increasing  w  in  (50),  '  becomes  small  (zero  for  the 

periodic  cases),  and  R  A  =:  1. 
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